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Abstract 

Let X n — {x^}^ =1 be a set of n points in the d-cube I d := [0, l] d , and $„ = {<ySj}" =1 a 
family of n functions in the space L q (I d ), < q < oo. We consider the approximate recovery 
in L q (I d ) of functions / on I d from the sampled values fix 1 ), f(x n ), by the linear sampling 
algorithm 

n 

L n (X n ,$ n J) := ^fix^ipp 
i=i 

Functions / to be recovered are from the unit ball of Besov type spaces of an anisotropic 
smoothness, in particular, spaces B^f of a "hybrid" of mixed smoothness a and isotropic 
smoothness (3. We constructed optimal linear sampling algorithms on special 

sparse grids X* and a family <£>* of linear combinations of integer or half integer translated 
dilations of tensor products of B-splines. We computed the asymptotic of the error of the 
optimal recovery. As consequences we obtained the asymptotic of optimal cubature formulas 
for numerical integration of functions from from the unit ball of these Besov type spaces. 

Keywords and Phrases Linear sampling algorithm • Cubature formula • Sparse grid • Besov 
type space of anisotropic smoothness • Quasi-interpolant representations by B-spine series. 

Mathematics Subject Classifications (2000) 41A15 • 41A05 • 41A25 • 41A58 • 41A63. 

1 Introduction 

The aim of the present paper is to investigate linear sampling algorithms and cubature formulas 
on sparse grids, and optimal sampling recovery and cubature of functions on the unit <i-cube 
I d := [0, l] d , having an anisotropic smoothness. 

Let X n = {x J }" =1 be a set of n points in I d , = {</?.,■ }™ =1 a family of n functions on I d and 
/ a function on I d . For approximately recovering / from the sampled values /(x 1 ), ...,f(x n ), we 
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define the linear sampling algorithm L n (X n , <£ n , •) by 

n 

L n (X n ,<S> n ,f) := X>(^>i- 
i=i 

Let -B be a quasi-normed space of functions on I d , equipped with the quasi-norm || ■ \\b- We measure 
the recovery error by ||/ — L n (X n , $ n , If W is a subset in B, to study optimality of linear 

sampling algorithms of the form for recovering / € W from n their values, we will use the 

quantity 

r n (W) B := inf sup ||/ - L n (X n , $ n , /)|| B . 

A general nonlinear sampling algorithm of recovery can be defined as 

R n (X n ,P n ,f) := P^/fcc 1 ), .»,/(*")), 

where P n : R n — > L q is a given mapping. To study optimal nonlinear sampling algorithms of 
recovery for / € W from n their values, we can use the quantity 

Qn(W) B := inf sup \\f - R n (X n , P n , f)\\ B - 

Let L q := L q (l ), < q < oo, denote the quasi-normed space of functions on I d , equipped with 
the gth integral quasi-norm |[ • || 9 for < q < oo, and the sup- norm || • for q = oo. We use the 
abbreviations: r n (W) q := r n (W) Lq and Q n (W) q := Q n (W) Lq - 

Recently, there has been increasing interest in solving approximation and numerical problems 
that involve functions depending on a large number d of variables. The computation time typi- 
cally grows exponentially in d, and the problems become intractable already for mild dimensions 
d without further assumptions. This is so called the curse of dimensionality. A classical model in 
attempt to overcome it which has been widely studied in literature is to impose certain anisotropic 
smoothness conditions on the function to be approximated and to employ sparse grids for con- 
struction of approximation algorithms. We refer the reader to [21 [131 (13 EH] f° r surveys and the 
references therein on various aspects of this direction. In the present paper, we focus our attention 
to sparse grids for problems of optimal sampling recovery and cubature. 

In sampling recovery and cubature problems for a function class W of anisotropic smoothness, 
the most important and challenging is the problem of constructing asymptotically optimal sampling 
algorithms L n (X n , •) and R n (X n , P n , •) in the sense of the quantities r n (W) b and g n (W) b and 
computing their asymptotic order. In solving this problem a key role play constructions of an 
appropriate set of sample points X n and of a series representation for functions having a given 
anisotropic smoothness. 

Sparse grids for sampling recovery and numerical integration were first considered by Smolyak 
|22j . For sampling recovery of periodic functions, he constructed the following grid of dyadic points 
in I d 

r(m) := {2- k s : k £ D(m), s G I d (k)}, 
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where D{m) := {k G Z d + : \k\! < m} and I d (k) := {s G 1 d + : < s { < 2 k \ i G [d]}. Here and in 
what follows, we use the notations: xy := (xiyi, XdUd)', 2 X := (2 Xl , 2 Xd ); \x\% := Yli=i \ x i\ f° r 
x, y G R rf ; [d] denotes the set of all natural numbers from 1 to d; Xi denotes the ith coordinate of 
x G R rf , i.e., x := (xi, It is remarkable that T(m) is a sparse grid of the size 2 m m d ~ 1 in 

comparing with the standard full grid of the size 2 dm . 

Temlyakov [231 [25] and Dinh Dung [8, 9j developed Smolyak's construction for studying the 
asymptotic order of r n (W) q for periodic Sobolev classes W® 1 and Nikol'skii classes Hp 1 having 
mixed smoothness a, as well their intersections, where 1 := (1, 1, 1) G R d . Recently, Sickel and 
Ullrich [20] have investigated r n (Upl) q for periodic Besov classes. For non-periodic functions of 
mixed smoothness 1/p < a < 2, this problem has been recently studied by Bungartz and Griebel 
[2] , Triebel |26] (d = 2) , Sickel and Ullrich |21] , using the mixed tensor product of piecewise linear 
B-splines and the grids of sample points T(m). Smolyak grids are a counterpart of hyperbolic 
crosses which are frequency domains of trigonometric polynomials widely used for approximations 
of functions with a bounded mixed smoothness. These hyperbolic cross trigonometric approxima- 
tions are initiated by Babenko [I]. For further surveys and references on the topic see O [25], the 
references given there, and the more recent contributions [20, 27 \. 

In the recent paper [12], we have studied the problem of sampling recovery of functions on I d 
from the non-periodic Besov class U^q, which is defined as the unit ball of the Besov space 
having mixed smoothness a. For various < p, 8,q < oo and 1/p < a < r, we proved upper 
bounds for r n (Upg) q which in some cases, coincide with the asymptotic order 

r n {U$) q x n-^/P-i/^+iog^-DV (1.2) 

where b = b(a,p,8,q) > and x+ := max(0, x) for x G R. By using a quasi-interpolant repre- 
sentation of functions / G g by mixed B-spline series, we constructed optimal linear sampling 
algorithms on Smolyak grids T(m) of the form 

L n (x*,$* n j) = E E /( 2 " fc j>fej» ( L3 ) 

k€D(m) jel d (k) 

where X* := T(m), $* := {ipk,j}keD(m),j£i d (k)' n '■= l r ( m )l x 2 m m d ~ 1 and ip k>j are explicitly 

(r) 

constructed as linear combinations of at most at most N B-splines Mf, ' for some N G N which is 

(r) 

independent of k,j,m and /, Mf, ' are tensor products of either integer or half integer translated 
dilations of the centered B-spline of order r > a. 

It is necessary to emphasize that any sampling algorithm on Smolyak grids always give the 
asymptotic of error as in (|1.2p with the logarithm term logg n, b > 0. Unfortunately, in the 
case when the dimension d is very large and the number n of samples is rather mild, the main 
term becomes log^ ^ & n which grows fast exponentially in d. To avoid this exponential grow 
we introduce other anisotropic smoothnesses and construct appropriate sparse grids for functions 
having these smoothnesses. Namely, we extend the above study to functions on I d from the classes 
Up f for a > 0,/3 G R, and U* g for a G R+, which are defined as the unit ball of the Besov type 

spaces and B^ e . The space B^ e and B^f are certain sets of functions with bounded mixed 
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modulus of smoothness. Both of them are generalizations in different ways of the space B^q of 
mixed smoothness a. The space B^ 9 is B^ \ for a = al. The space B"f is a "hybrid" of the space 
B^\ and the classical isotropic Besov space £?£ R of smoothness B. The parameter a governs the 
mixed smoothness, whereas 8 in B® g governs the isotropic smoothness. It coincides with B^ for 
B = 0, and with B^ for a = 0. 

The space B p ' e is a Besov type generalization of the Sobolev type space H a ^. The latter space 
has been introduced in |15tll6| for solutions of elliptic variational problems in high dimensional set- 
tings. By use of tensor-product biorthogonal wavelet bases, the authors of these papers constructed 
so-called optimized sparse grid subspaces for finite element approximations of the solution having 
H a -,P '-regularity, whereas the approximation error is measured in the norm of classical isotropic 
Sobolev space i? 7 . 

In the present paper, with some reasonable restrictions, the asymptotic orders of r n (Upf) q , 

Qn{Up'0) q for the case 8 ^ 0, and r n {U^ e ) q , Q n {Up d ) q for the case ai < d2 < ... < are completely 
computed. In particular, if < p, 6, q < oo, we prove for a + 8 > 1/p and 8 < 0, 

rn(U;f) q x en (U a J) q * n -~-W/p-i/<0 +> (1.4) 

and for a\ > 1/p, 

r n (U a pfi ) q x g n (U; fi ) q x n -x+(l/P-l/*) + . (1.5) 

It is remarkable that these asymptotic orders do not contain any exponent in d and moreover, do 
not depend on d. For a set A C Z^, we define the grid G(A) of points in I d by 

G(A) := {2- k s : k € A, s € I d (k)}. (1.6) 

For the quantities of optimal recovery in (jl.4p and (II. 5p . asymptotically optimal linear sampling 
algorithms of the form 

L n (X*,$* n ,f) = E f( 2 ~ k j)^j (1-7) 

fceA jei d (k) 

are constructed where X* := G(A), $* := {"^fejIfceA, j&i d (k) an d ^fcj are the same as in (11.31) . 
The set A which is parameterized as A(£) or A'(£) in ^ > 0, is specially constructed for each 
class of Up f and U® 9 , depending on the relationship between < p, 9, q < oo and a, 8 or between 
< p, 9,q < oo and a, respectively. The grids G(A(£)) or G(A'(£)) are sparse and have much 
smaller number of sample points than the standard full grids, and give the same error of the 
sampling recovery on the latter ones. The asymptotically optimal linear sampling algorithms 
L n (X*, <]>* , ■) are based on quasi-interpolant representations by B-spline series of functions in 
spaces Bpf and B^q. We also compute the asymptotic order of f n {Upf) B i T and g n (Up'g) B J T for 
the case B 7, and construct corresponding asymptotically optimal linear sampling algorithms of 
the form ()1.T|) f)1.6[) . As consequences of these results, we obtain asymptotically optimal cubature 
formulas for functions from these classes and asymptotic order of the error of the approximate 
integrations by them. 

We are restricted to compute the asymptotic order of r n (Upf) q and Q n (Upf) q with respect 
only to n when n — > 00, not analyzing the dependence on d. Recently, in [13] Kolmogorov n-widths 
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d n (U, H 1 ) and e-dimensions n £ (U, H 1 ) in space H" 1 of periodic multivariate function classes U have 
been investigated in high-dimensional settings, where U is the unit ball in H a, ° or its subsets. We 
computed the accurate dependence of d n (U,H' y ) and n^U^H 1 ) as a function of two variables n, 
d or e, d. Although n is the main parameter in the study of convergence rate with respect to n 
when n —> oo, the parameter d may affect this rate when d is large. It is interesting and important 
to investigate r n (U, H 1 ), g n {U,H rY ) or more generally r n (U"'g , Bq >T ), Q n {U^f , Bq )T ) and cubature 
in these high-dimensional settings. We will discuss this problem in a forthcoming paper. 

The present paper is organized as follows. In Section we give definitions of Besov type 
spaces of of functions with bounded mixed modulus of smoothness, in particular, spaces 

Bp'g and B^ e , and prove theorems on quasi-interpolant representation by B-spline series, and 
relevant discrete equivalent quasi-norms for these spaces. In Section[3j we construct linear sampling 
algorithms on sparse grids of the form (II. 7ft for function classes U^f and U® e , and prove upper 
bounds for the error of recovery by these algorithms. In Section [H we prove the sparsity and 
asymptotic optimality of the linear sampling algorithms constructed in Section [3l for the quantities 
f n (Up'g) q , ^n(Up d ) q , Qn{U^f) q , Qn{Up g) q and establish their asymptotic orders. In section [5j we 

extend the investigations of Sections [3] and H] to the quantities T n (U^f) B i T and Q n (Upf)Bj T - 
In Section [6l we discuss the problem of optimal cubature formulas for numerical integration of 
functions from U^f and U^ e . 



2 Function spaces and B-spline quasi-interpolant representations 

Let us first introduce spaces B^ g of functions with bounded mixed modulus of smoothness and 

Besov type spaces B^f , B^ of functions with anisotropic smoothness and give necessary knowl- 
edge of them. 

Let G be a domain in R. For univariate functions / on G the rth difference operator A r h is 
defined by 

AUf,x):= j2(-iy^( r )f(x + jh). 

3=0 W 

If e is any subset of [d], for multivariate functions on G d the mixed (r, e)th difference operator A^' e 
is defined by 

where the univariate operator is applied to the univariate function / by considering / as a 
function of variable Xi with the other variables held fixed. For a domain G in R d , denote by L p (G) 
the quasi-normed space of functions on G with the pth integral quasi-norm || • \\ p q for < p < oo, 
and the sup norm || ■ \\oo,G for p = oo. Let 

^r(/^W : = SU P \\ A h e (f)\\ P ,G d (h,e), * G G<i > 
\hi\<ti,i£e 
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be the mixed (r, e)th modulus of smoothness of /, where G d (h, e) := {x G G d : X{,Xi + rhi GG, i£ 
e} (in particular, wf (f,t) p>G d = ||/|| P)G d). 

For x,x' € M. d , the inequality x' < x (V < x) means x', t < Xi (x[ < Xj), i G [d]. Denote: 
M + := {x G E : x > 0}. Let Q : — >• M + be a function satisfying conditions 

n(t) > 0, f > 0, t G Rj, (2.1) 

n(t) < CO(t'), t<t', M'eRt (2.2) 
For a fixed 7 G R+, 7 > 1, there is a constant C" = C'i^f) such that for every A G R+ with A < 7, 

n(\t) < c'n(t), teR d + . (2.3) 

For e C [d], we define the function f2 e : R^ — > R+ by 

fi e (t) := 0(i e ), 

where t e G R+ is given by tj = tj if j G e, and i| = 1 otherwise. 



If < p, 9 < 00, we introduce the quasi-semi-norm \ f\ B n i e \ for functions / G L p (G d ) by 



f K(/,t) P)G ,/H e (t)f [Tf: 1 ^ , 0<oo, 
7ld / (2.4) 



SUp Wr(/,*)p,G d / fi e(*), 61 = OO, 



(in particular, |/| B n 



Another alternative definition of the quasi-semi-norm \f\B n e ( e ) * s obtained by replacing the 

integral or supremum over I d in ([MD by one over Mi_. In what follows, we preliminarily assume 
that the function £1 satisfies the conditions (|2.1|) - (|2.3|) . 

For < p, 9 < 00, the Besov type space B^ e (G d ) is defined as the set of functions / G L p (G d ) 
for which the quasi-norm 

eC[d] 

is finite. Since in the present paper we consider only functions defined on I d , for simplicity we drop 
the symbol I d in all the above notations. 

Lemma 2.1 Let < p, 9 < 00. Then there holds true the following quasi-norm equivalence 

ii/Hbo, - Mf) ■= E ( E {<{f^- k ) P m- k )} ) , 

with the corresponding change to sup when 9 = 00. 
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Proof. This lemma follows from properties of mixed modulus of smoothness u^(f,t) p and the 
properties (|2.ip ~ (|2.3p of the function CI. We prove it for completeness. The lemma will be proven 
if we show that for every e C [d] , 



with the corresponding change to sup when 6 = oo. Let us prove this semi-norms equivalence for 
instance, for e = [d], 1 < p < oo and < 9 < oo. The general case can be proven in a similar way 
with a slight modification. 

Put D(k) := {x e R'l : k < x < k + 1} and use the abbreviation u r (f, -) p := uif\f, -) p . By 
(|2HJ)-(|23D we have 

Q(2~ x ) x ft(2" fc ), x G D(jfe), A; G Z^.. (2.5) 
From the monotonicity of u r (f, •) in each variable and the inequality 

d 

Wr(/,ci) p < [J(l + Cj -) r a; r (/,t) p , crtf, c> 0, 

we obtain 

^(7,2-*), x ^ r (/,2- fc ) p , xeD(k), fcezf. (2.6) 
Setting /(A;) := {tel d : 2-^ < t < 2~ k }, by ||23J) and (USD we have 



i^ 9 ([«a) = e/ K(/,tvn(t)}*n^ 

V/ K(/,2-) p /0(2- 2 )}^ x ^K(/,2~ i yo(2- i )} ( 



Let us define the Besov type spaces B^f and B£ e of functions with anisotropic smoothness as 



particular cases of B^ 

Let a G R+ and j3 G R with a + /3 > 0. We define the space B^'g as follows. 



VP 



B$ := B£„ where «(t) = < 



i al inf tf , /3 > 0, 

(2.7) 

i al sup t£, /3 < 0. 
ie[d] 



The definition (|2.7p seems different for /3 > and /3 < 0. However, it can be well interpreted in 
terms of the equivalent discrete quasi-norm B\{f) in Lemma 12. 11 Indeed, the function Cl in (12. 7p 
for both (3 > and /3 < satisfies the assumptions ()2. 1 1) f)2.3[) and moreover, 



1/0(2""*) = 2 a ^ 1+ PW°°, x£^ d 



+ • 
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where \x\oo := max^gui \xj\ for (see the proof of Theorem 12. 2p . Hence, by Lemma [2.11 there 

holds true the following quasi-norm equivalence 

Wf\\B«J ~ E ( E {2 a|fc|l+/3|fc| -<(/,2- fc ) p } J (2.8) 

with the corresponding change to sup when 6 = oo. The notation 5^ ^ becomes explicitly reason- 
able if we take the right side of (|2.8j) as a definition of the quasi-norm of the space B^f . 

The definition of B^f includes the well known classical isotropic Besov space and its mixed 
smoothness modifications. Thus, we have B^f = B^\ for j3 = 0, and B"'g = B^ e for a = 0, where 
Bp is the classical isotropic Besov space of smoothness (3. 

For a G M+, we define the space B^ 9 of mixed smoothness a as follows. 

B a vfi := Bp e , where O(i) = t a , t G . (2.9) 

Next, we introduce quasi-interpolant operators for functions on For a given natural number 
r, let M be the centered B-spline of order r with support [— r/2,r/2] and knots at the points 
— r/2, — r/2 + 1, r/2 — 1, r/2. Let A = {A(s)} :) ' g p( At ) be a given finite even sequence, i.e., A(— j) = 
X(j), where -P(/i) := {j G Z : |j| < /i} and fi > r/2 — 1. We define the linear operator Q for 
functions / on M by 

Q(/,x) := J2 A (f,s)M(x-s), (2.10) 

where 

A(/,*):= ^ A(i)/( S -j). (2.11) 

The operator Q is local and bounded in C(M) (see [3l p. 100-109]), where C(G) denotes the 
normed space of bounded continuous functions on G with sup-norm || • ||c(G)- Moreover, 

||0(/)||CCR) < l|A||||/||ocR) 

for each / G C(R), where ||A|| = Ej £ p( M ) An operator Q of the form (l2TTUj) - (f2TTTD 

reproducing V r — i, is called a quasi-interpolant in C(R). 

If Q is a quasi-interpolant of the form (|2.10p - (|2.1ip . for h > and a function / on M, we define 
the operator Q(-;h) by 

Q(/;/i) := cr h oQ o a 1/h (f), 
where ah{f,x) = f(x/h). From the definition it is easy to see that 

Q(f,x;h)= Y J Hf,k;h)M(h- 1 x-k), 
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where 



A(f,k;h):= Yl Kj)f(h(k-j)). 



The operator Q(-;h) has the same properties as Q: it is a local bounded linear operator in 
C(R) and reproduces the polynomials from V r —x- Moreover, it gives a good approximation for 
smooth functions [4, p. 63-65]. We will also call it a quasi-interpolant for C(K). However, the 
quasi- interpolant Q(-;h) is not defined for a function / on I, and therefore, not appropriate for an 
approximate sampling recovery of / from its sampled values at points in I. An approach to con- 
struct a quasi-interpolant for functions on I is to extend it by interpolation Lagrange polynomials. 
This approach has been proposed in [10] for the univariate case. Let us recall it. 

For a non-negative integer k, we put Xj = j2~ k ,j £ Z. If / is a function on I, let U k (f) and 
V k (f) be the (r-l)th La grange polynomials interpolating f at the r left end points Xo, xi, x r _i, 
and r right end points x 2 k_ r+1 , x 2 k_ r+3 , ...,x 2 fc, of the interval I, respectively. The function f k is 
defined as an extension of / on M by the formula 



fkix) :-- 



U k (f,x), x<0, 
f(x), 0<x<l, 

y k (f, x ), x >i. 



If / is continuous on I, then f k is a continuous function on M. Let Q be a quasi-interpolant of the 
form dZTOD - (gJT} in C(R). Put Z + := {k € Z : k > -1}. If k G Z+, we introduce the operator 
Qfc by 

Qk(f,x):= Q(/ fc ,x;2- fc ), and Q_i(/, x) := 0, s 6 I, 



for a function / on I. We have for k G 



Qk(f,x) = ^2 a kjS (f)M kjS (x), Vx G I, 

where ./(A;) := {s G Z : — r/2 < s < 2 k + r/2} is the set of s for which M k s do not vanish 
identically on I, and the coefficient functional a k<s is defined by 

a M (/):= A(/ fc , S ;2- fc )= ^ A(j)/ fc (2- fe ( S - j)). 

lil<A« 



Put := {£; G Z+ : k{ > —1, i G [d]}. For A; G Z^, let the mixed operator Q k be defined by 

d 

Qk-.= Y[Qki, (2-12) 

1=1 

where the univariate operator Q^. is applied to the univariate function / by considering / as a 
function of variable Xj with the other variables held fixed. 
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We define the integer translated dilation Mf- S of M by 

M k<s (x) := M{2 k x - s), G Z+, sGZ, 
and the d- variable B-spline M ktS by 

d 

M k A x ) ■= YlM kijSi ( Xi ), keZ d + , se Z d , (2.13) 
1=1 

where Z + is the set of all non-negative integers, Z^ := {s G % d : Sj > 0, i G [d]}. Further, we 
define the half integer translated dilation of M by 

M^ s {x) ■= M(2 k x - a/2), fc G Z+, s G Z, 
and the d- variable B-spline M£ s by 

d 



M fc %(x) := ][[ M fcW^)> fc G Z^., s G Z d . 



i=l 

(r) 

In what follows, the B-spline M will be fixed. We will denote Mt ' := M kjS if the order r of M is 

(r) 

even, and s := M£ g if the order r of M is odd. 
We have 

Qk(f,x) = Yl "kAfWkM* VxGi d , 

seJ d {k) 

where M M is the B-spline defined in (|2TT3|) . J d (fc) := {s G Z d : -r/2 < < 2 k ' + r/2, i G [d]} 
is the set of s for which Mfc s do not vanish identically on I , 

a>k,s(f) '■= a ki,si(a>k2,s 2 (- a k d ,s d (f)))> ( 2 - 14 ) 

and the univariate coefficient functional a kiS . is applied to the univariate function / by considering 
/ as a function of variable x% with the other variables held fixed. 

The operator Q k is a local bounded linear mapping in C(I d ) for r > 2 and in for r = 1, and 
reproducing V d _\ the space of polynomials of order at most r — 1 in each variable ccj. In particular, 
we have for every / G C(I d ), 

llft(/)l|oo<C||A|| d ||/|| c(Id) . (2.15) 

For k G Z+, we write k — > oo if k{ — > oo for i G [d]). 
Lemma 2.2 We have for every f G C(I d ), 

\\f-Qk(f)\\oo < C Yl ^r(/> 2-fe )oo, (2.16) 

eG[d], e^0 

and, consequently, 

||/-Qfc(/)IU^0, fc^oc. (2.17) 
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Proof. For d = 1, the inequality (|2.16p is of the form 

||/-Qfc(/)||oc < Cu r {f, 2- fc )oo. (2.18) 

This inequality is derived from the inequalities (2.29)-(2.31) in [11] and the inequality (|2.15p . For 
simplicity, let us prove the the inequality (|2.16|) for d = 2 and r > 2. The general case can be 
proven in a similar way. Let / be the identity operator and k = (k\, k%). From the equation 

I-Qk = (I - Qfci) + {i - Qk 2 ) - {I - Qki)(l - Qk 2 ) 

and the inequality ()2. 18|) applied to / as an univariate in each variable, we obtain 

11/ - QuU)U < W- QkMWoo + - Qfe,)(/)||oo + ll(i - Q kl )(i - Qk 2 )(f)\\oo 
« cjP(f, 2-^)00 + u4 2 >(/, 2-^)00 + 4 2 i(/, 2- fc ) 00 . 



If r is a number such that < r < mm(p, 1), then for any sequence of functions {gt} there is 
the inequality 

|E^r ^ En^iip- ( 2 - 19 ) 



Let J r d (/c) := J d {k) if r is even, and J d (k) := {s e Z d : -r < Si < 2 fcl+1 + r,i € [d]} if r is odd. 
Notice that J d {k) is the set of s for which M^ r j do not vanish identically on I d . Denote by T*f(k) 

the span of the B-splines M> r j, s £ J d (k). If < p < 00, for all k 6 and all g € S^(A;) such 
that 

9= E a * M K' ( 2 - 2 °) 

there is the quasi-norm equivalence 

\\g\\ p X 2-\ k ^\\{a s }\\ Ptk , (2.21) 

where 

IIKIIU* := ( E l^l P ) 7 



with the corresponding change when p = 00. 

Let the operator q^, k € Z+, be defined in the manner of the definition (I2.12P by 

d 



We have 



(Ik := H(Q ki ~ Q ki -l) ■ (2.22) 

i=l 

Qk = E^'- ( 2 - 23 ) 
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From (|2.23j) and (|2. 17j) it is easy to see that a continuous function / has the decomposition 

/ = E %(/) 

with the convergence in the norm of L^. 

From the definition of (|2.22p and the refinement equation for the B-spline M, we can represent 
the component functions %(/) as 

**(/)= E c £(/)<i> (2-24) 
where eft are certain coefficient functionals of /, which are defined as follows (see [12] for details) . 



We first define c[ r ^ for univariate functions (d = 1). If the order r of the B-spline M is even, 



<£(/) 



: = «m(/)-4, s (/), ^>0, (2.25) 
where 

a' M (/) := 2~ r+1 Y, (%k-i,m(f), k > 0, <.(/):= 0. 
(mj)eCV(M) 

and 

C r (fc,s) := {(m,j) : 2m + j - r/2 = s, m G - 1), <j < r}, fc > 0, C P (0,s) := {0}. 

If the order r of the B-spline M is odd, 

*0, fc = 0, 

a kjS/2 (f), k > 0, s even, 

l2" r+1 EMeftM G)»*-i,m(/). * > °> s odd > 
where 

C r (k,s) := {(m, j) : 4m + 2j - r = s, m G J(fc - 1), < j < r}, k > 0, C r (0,s) := {0}. 

In the multivariate case, the representation (|2.24p holds true with the eft which are defined in 
the manner of the definition of (|2,14p by 

$(/) = ctiML(--4ism- (2-26) 

Thus, we have proven the following 
Lemma 2.3 Every continuous function f on I d is represented as B-spline series 

f = E »(/) = E E %MWi% (2.27) 

converging in the norm of L^, where the coefficient functionals cft s (f) are explicitly constructed 
by formula (|2.25|) - (|2.26p as linear combinations of at most N function values of f for some N G N 
which is independent of k,s and f. 
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We now prove theorems on quasi- interpolant representation of functions from B^ d and B^'g , 
Bp 6 by series (|2.27p satisfying a discrete equivalent quasi-norm. We need some auxiliary lemmas. 

Let us use the notations: x + := ((xi)+, (xd)+) for x G M. d ; := {s E : Sj = 0, j ^ e} 

and N d (e) := {s G : Sj > 0, i£e, Sj = 0, i ^ e} for e C [d] (in particular, N d (0) = {0} and 
N d ([d]) = N d ). We have N d (u) n N d (?;) = if « / v, and the following decomposition of Z^_: 

= |J N d (e). 

eC[<i] 

Lemma 2.4 ([12]) Let < p < oo and r < min(p, 1). Then for any f G C{l d ) and k G N d (e), 
i/iere /io/ds £/ie inequality 

( \ 1 / T 

\Mf)\\ P < C£ £ {2l- fc I^V(/,2- 

^^Z^O), s>k 

with some constant C depending at most on r,fj,,p,d and ||A||, whenever the sum in the right-hand 
side is finite. 

Lemma 2.5 Let < p < oo, < r < min(p, 1), 5 = min(r, r — 1 + Xei g G L p be represented 
by the series 

9 = E 9k, £r(*0 
converging in the norm of L^. T/ien /or any /c G Z+(e), i/iere holds the inequality 

/ \ l / T 

u e r (g,2- k ) p < C 



e {^ i{k ~ s)+ H 9 s\\ P y) 



wrai/i some constant C depending at most on r,fi,p,d and ||A||, whenever the sum on the right-hand 
side is finite. 

Proof. This lemma can be proven in a way similar to the proof of [12\ Lemma 2.3]. □ 

Let < p, 6 < oo and ip : Z+ — > R. If {gk}kez^ is a sequence whose component functions gk 



l/e 



are in L p , we define the "quasi-norm" ||{gfc}|| b </' by 

Winn,,*, ■■= (e { 2m \\9k\\ P y) 

kez^_ 

with the usual change to a supremum when 9 = oo. When {gk}kez d * s a positive sequence, we 
replace ||ofc|L by \gk\ and denote the corresponding quasi-norm by ||{<7fc}||,j/>- 

We will need the following generalized discrete Hardy inequality (see, e.g, [5] for the univariate 
case with tp(k) = ak, a > 0). 
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Lemma 2.6 Let {ak\k&'E d an ^ {bk}k£Z d be ^ wo positive sequences and let for some M > 0, r > 
0, 5 > 

( r\ 1A 

bk < M Yl (2 <5|(s_fc)+|l a s ) T . (2.28) 
^sez£ J 

Let the function ip : — > R satisfy the following. There are numbers q,C2 6 M, e > and 
< £ < 5 such that 

ip(k)-e\k\i < ip(k') -e\kf\i + a, k < k', k,k' G Z+, (2.29) 

and 

^(fc)-Cl*|i > ^00 - - 02, k <k', k,k' GZf, (2.30) 

TTien for < 9 < oo, there holds true the inequality 

\\{b k }\\ bt < CM\\{a k }\\ b * (2.31) 

with C = C(ci, C2, e, (5, 0, d) > 0. 

Proof. Because the right side of (12.28P becomes larger when r becomes smaller, we can assume 
t < 9. From (|2.28j) we have 

b k < M ^ B fc (e), fceZf, (2.32) 

eC[d] 

where 



B*(e):= 2^)1^ ^ (2 'a„ 



5<5|s(e)|i. 

s&Z(e,k) 



T 



and 

Z(e, fc) := {s £ Z{ : Sj < fcj, j G e; Sj > fcj, j e}. 

For e C [d] and s G Z d , let e := [d] \ e and s(e) G Z d be defined by s(e)j = Sj if j G e 
and s(e)j = if j e. Take numbers with the conditions < e' < e, C < C < $ an d 

t/0 + t/9' = 1, respectively. Applying Holder's inequality with exponents 9/t,9'/t, we obtain 

/ \ V» / \ I/"' 

Sfc(e) < 2~ 5|fc(e)|1 ( ^ (V' |s(e)|l+e ' |s(5)|l a s M ( ^ ^-Ol^li-c'We)!/ 

\seZ(e,fe) / \se^(e,fc) 



K s£Z(e,k) 



< 2 — C[fcCe)U-*'|fc(e)li j Y (V' |s(e)|l+e ' |s(g)|l a. 

\ sGZ(e,fc) 



1/0 

6» 
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Hence, 



seZ(e,k) 

(2 33) 

< ^ 2 9K ' |s(e)|l+e ' |s(s)|l) af V 2 e(,/,(fe) - c ' |A;(e)|l - e ' |fc(5)|l) , 

s&,\ keX(e,s) 

where 

X(e, s) := {k € : fcj > Sj, j € e; fcj < Sj, j e}. 

By (gSSD and (jZ5D|| we have for fc G X(e, s), 

^(fc) = ^(fc)-Cl*|l + C(|fc(e)|i + |*:(e)|i) 

< ^(s(e), fe(e)) - C(|«(e)|i + |fc(e)|i) + C(|A:(e)|i + |fc(e)|i) 
= ^( s ( e ),A:(e))-C| S (e)|i + C|A : (e)|i, 

and 

^(s(e),k(e)) = iP(s(e),k(e)) - e(|s(e)|i + [fc(e)[i) + e(| S (e)|i + [fc(e)|i) 
< ^( S (e), s(e)) - e(| S ( e )|! + | S (e)|i) + e([s(e)[i + |fc(e)|i) 
= ^0) - e|s(e)|i + e|fe(e)|i. 

Consequently, 

^(fc) - C'|fc(e)|i - e'|fc(e)|i < Ms) - C|«(e)|i - e|a(e)|i - (C' - C)l*(e)|i + (e - e')|fc(e)li> 
and therefore, we can continue the estimation (|2.33p as 

sez d + kex(e,s) 

< 2 e W s )+(C-0\s(e)\i-(e-e')\s(e)\ la e 2 0(-(C'-C)| S (e)|l+(e-e')K§)|l 

= £ 2*««f = ||{a*}||;,. 

Ug 

Hence, by (^321) we prove (|2^B . □ 

We now are able to prove quasi-interpolant B-spline representation theorems for functions from 
and B^'g, B^g. For functions / on I d , we introduce the following quasi- norms: 

1/6 



B*(f) ■= (E {!k fe (/)ll P /^(2- fc )} 

fla(/) := ( E {2- |fe|l/p ||{c2(/)}|| M /0(2- fe )} 
Observe that by (|2.2ip the quasi-norms B%(f) and B^(f) are equivalent. 
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Theorem 2.1 Let < p, 9 < oo and £1 satisfy the additional conditions: there are numbers 
p, p > and Ci, C2 > such that 

n(t)t-'" 1 < CiO(0£'" Ml , t<t', t,t' el d , (2.34) 

n(t)r pl > c 2 n(t')t'- pl , t>t',t,t'£i d . (2.35) 

Then there hold the following assertions. 

(i) If p > 1/p and p < r, then a function f £ B^ s can be represented by the B-spline series 
(|2,27p satisfying the convergence condition 

B2U) « \\f\\ B n. (2.36) 

P ,e 

(ii) If p < min(r, r — 1 + 1/p), then a continuous function g on I d represented by a series 

9 = E»=E E c ^Kl 



satisfying the condition 



B ^9):= (E {\\gk\\ P M2- k )yy 9 < 



OO, 



belongs the space B^ e . Moreover, 

\\g\\ B n < B 4 {g). 

(Hi) If p > 1/p and p < min(r, r — 1 + 1/p), then a continuous function f on I d belongs to the 
space BpQ if and only if f can be represented by the series (|2.27p satisfying the convergence 

condition (|2.36p . Moreover, the quasi-norm is equivalent to the quasi-norm B 2 (f). 

P ,g 

Proof. Put 4>(x) := log 2 [l/fi(2 -a: )]. Due to ^Mj) - <^35h . the function <f> satisfies the following 
conditions 

<j)(x) — p\x\\ < 4>{x') — p\x'\i + log 2 Ci, x<x', x,a;'elj| (2.37) 

and 

<j){x)-p\x\i > 4>{x') - p\x'\ x - log 2 C 2 , x < x', x,x' e R%. (2.38) 

We also have 

B ^f) = E ( E {^(2 fc K(/,2- fc ) p )} 9 ) , (2.39) 
with the corresponding change to sup when = oo. Fix a number < r < mm(p, 1). 
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Assertion (i): From ([2.37P we derive fi\k\% < 4>(k) + c, k G Z+, for some constant c. Hence, by 
Lemma (|2.ip and (|2.39j) we have 

\\f\\ B »i < C\\f\\ B a f G BpQ, 

for some constant C. Since for \i > 1/p, IJp 1 is compactly embedded into C(I d ), by the last 
inequality so is B^ g . Take an arbitrary / G -B^V Then / can be treated as an element in C(I d ). 
By Lemma [2,31 / is represented as B-spline series (|2.27p converging in the norm of L^. For k G Z^_, 
put 



b k :=2\ k ^\\ qk (f)\\ p , a k := ( £ {^^(f, 2" fc ) p } T ) 



1/T 



if G N d (e). By Lemma |2H we have for k G Zi, 





l/T / 




i » 







l/r 



for a fixed 5 > p + Let the function ^ be defined by ip(k) = <ft(k) — \k\i/p, k G lA. By the 
inequality \i > 1/p, (|2.37p and (|2.38[) . it is easy to see that 



^(fc)-e|fc|i < W) -e\k'\i +log 2 Ci, k<k',k,k'e 



J +' 



and 



i/)(k)-C\k\i > ip(k') — CW\i — log 2 C 2 , k<k', k,k' eZ%, 
for e < p — 1/p and Q = p + 1/p. Hence, applying Lemma [2761 gives 

B*tf) = \\{b k }\\ h * < C\\{a k }\\ b * ~ Bi(f) x . 



Assertion (ii): For k G Zi, define 



6fc :- 



1A 



Ofc := 



.■hp 



if A; G N d (e). By Lemma [231 we have for any A; G Z+(e), 

^( 5 ,2- fc ) p < C 3 \J2 {^ Slik - s)+ll \\9s\\ P 



l/r 



where 5 = min(r, r — 1 + 1/p). Therefore, 

b k < C 4 ^ (2 <5 K fc - s )+l 1 a s ) | , k G 
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Taking ( = p and < e < n, we obtain by ([2T37|) and ([2738]) 

(t){k)-e\k\i < 4>{k') -e|A/|i + log 2 Ci, k < k', k,k' e 

and 

0(fc)-C|fc|i > 0(^0 -Cl^'li -log 2 C 2 , k<k', k,k' 
Applying Lemma 12.61 we get 

||<7|U X B x (g) x ||{6 fc }|| , < C\\{a k }\\ b , = B 4 (g). 
Assertion (ii) is proven. 

Assertion (iii): This assertion follows from Assertions (i) and (ii). □ 
From Assertion (ii) in Theorem 12.11 we obtain 

Corollary 2.1 Let < p, 9 < oo and Q satisfy the assumptions of Assertion (ii) in Theorem \2.1\ 
Then for every k £ 

\\g\\ B s\ « \\g\\ P /^- k ), g € E d (k). 
Theorem 2.2 Ze£ < p, < oo and a G R+, {3 € R suc/i f/iaf 

l/p < min(a, a + fi) < max(a, a + 0) < r. 
Then there hold the following assertions. 

(i) A function f € B^'fi can be represented by the mixed B-spline series (|2.27p satisfying the 



convergence condition 



B2(f) 



( 



£ {^ lkll+mac \\Qk(f)\\ P } d \ « \\f\\ Kf - (2.40) 



£2 



(m,) if in addition, max(a,a + /3) < min(r, r — 1 + l/p), i/ien a continuous function f on I d 
belongs to the space B^f if and only if f can be represented by the series (|2.27p satisfying 
the convergence condition ()2.40p . Moreover, the quasi-norm ||/|| Rt *„s is equivalent to the 

quasi-norms i?2(/)- 

Proof. As mentioned above, for 9. as in ([121), we have l/tt{2~ x ) = 2 a l I 'l 1+/3 l a; l° , x G R^_. By 
Theorem 12. 1[ the assertion (i) of the theorem is proven if the conditions (|2.ip ~ (|2.3p and (I2.34p - 
(|2.35p with some p > l/p and p < r, are verified. The condition (|2.ip is obvious. Put 

(t>(x) :=log 2 {l/n(2- x )} = a|s|i+ 014,0, x eR d + . 

Then the conditions ()2.2j) — (|2.3|) and ([2.34p f[2.35|) are equivalent to the following conditions for the 
function <f>, 

<f>(x) < <f)(x')+log 2 C, x <x', x,x' £R d + ; (2.41) 
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for every b < log 2 7 := (log 2 71 , log 2 7^), 

<p(x + b) < <j){x) + log 2 C, + (2.42) 

<j>{x) -fi\x\i < 4>{x') - p\x |i + log 2 Ci, x < x', x, x' G M^; (2.43) 
0(x)-p|x|i > 0(a/) -p|x'|i -log 2 C 2 , x < x', x, x' G (2.44) 

We first consider the case /3 > 0. Take and p with the conditions 1/p < p < a and p = a + /3. 
The conditions (|2.4ip ~ (|2.42p can be easily verified. From the inequality a — p > and the equation 

(f>(x) - (i\x\i = (a — [i)\x\i + /3|x|oo, x G R+. ( 2 -45) 

follows ([2331) • We have 

0(;c) — p\x\i = /3(|x|oo — |x|i) = — /3 min \^ Xj, x G Ml. 

Hence, we deduce (|2.44|) . 

Let us next consider the case f3 < 0. The condition (|2.42h is obvious. Take p and /? with the 
conditions 1/p < p < a + (3 and a < p < r. Let x < x', x,x' E R+. Assume that |x|oo = Xj and 
l^'loo = x'j,. By using the inequalities xy < Xj < x'j < x'j, and a — p>a + /3 — p>0 from (12.45P 
we get 

cf)(x) — p\x\i = (a — p) X{ + (a + (3 — p)xj + (a — p)xj> 
< (a — p) x\ + (a + /3 — /i)x^-/ + (a - /u)x^ 
= </>(x') — /i|x'|i. 

The inequality (|2.43p is proven. The inequality (|2,4ip and (|2.44p can be proven analogously. 
Instead the inequalities a — p > a + /3 — p > 0, in the proof we should use a > a + (3 > and 
a + j3 — p < a — p < 0, respectively. Thus, the assertion (i) is proven. 

The assertion (ii) can be proven in a similar way. □ 

Theorem 2.3 Let < p, 8 < 00 and 7 G M+ such that 

d/p < 7 < r. 

Then there hold the following assertions. 

(i) A function f G Bj can be represented by the mixed B-spline series (|2.27p satisfying the 
convergence condition 

( 

B 2 (f) = 



£ {2^\Mf)\\p} \ « ll%V (2-46) 
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(ii) If in addition, 7 < min(r, r — 1 + 1/p), then a continuous function f on I belongs to the 
space Bq <T if and only if f can be represented by the series (|2.27p satisfying the convergence 
condition (|2,46p . Moreover, the quasi-norm \\f\\gy is equivalent to the quasi-norms I?2(/)- 



Proof. This theorem can be proven in way similar to the proof of Theorem 12.21 with a slight 
modification. In particular, in the proof of Assertion (i), the condition 1/p < min(a, a + j3) is 
replaced with the condition d/p < (3. □ 

Theorem 2.4 Let < p, 6 < 00 and a G such that 

1/p < mina,- < maxa,- < r. 
je[d] je[<q 

Then there hold the following assertions. 



(i) A function f £ B®g can be represented by the mixed B-spline series (|2.27p satisfying the 
convergence condition 



B2V) 



£ {2^ fc )||g fc (/)|| p } 6 I « Hjiifl., 



(2.47) 



(ii) If in addition, maxj g M dj < min(r, r — 1 + 1/p), then a continuous function f onI d belongs to 
the space B^ d if and only if f can be represented by the series (|2.27p satisfying the convergence 
condition (|2.47p . Moreover, the quasi-norm is equivalent to the quasi-norms B2(f). 



Proof. For Q as in (|2,9p . we have 

l/n(2~ x ) = 2^ x \ x G m d + . 

One can directly verify the conditions (|2.ip - (|2.3p and the conditions (|2.34p - (|2.35p with 1/p < fj, < 
m i n je[d] a j an d p = maxj g [^ aj, for Q defined in ()2.9p . Applying Theorem I2.1f i). we obtain the 
assertion (i). 

The assertion (ii) can be proven in a similar way. □ 

Remark Theorem 12.41 for a = al and Theorem 12.21 for (3 = coincide. This particular case has 
been proven in |12j . 



3 Sampling recovery 

Let AcZf be given. We define the operator for functions / on I d by 

RA(f):= !>(/) = E E c 2(/)<i, (3-1) 
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and the grid G(A) of points in I d by 

G(A) := {2- k s : k G A, s G 
Denote by M r d (A) the set of B-spines M^j, fe G A, a G J r d (fc). 

Lemma 3.1 The operator defines a linear sampling algorithm of the form (jl.ip on the grid 
G(A). More precisely, 

keA jei d (k) 

where X n := G(A), 3> n := {V^jlfcgA, je/ d (fc)> 



7? 



:=|G(A)| = ^H( 2% + 1 ) 

ieAj=l 



and ^ are explicitly constructed as linear combinations of at most N B-splines M^ r j G M d (A) 
/or some JVgN which is independent of k,j, A and f . 

Proof. This lemma can be proven in a way similar to the proof of [121 Lemma 3.1]. □ 



Let i/j : 1% -> R 
quasi-norm is finite 



. Denote by the space of all functions / on I d for which the following 



B 



{^ {k) \\Qk(f)\\ P } e 



k& d + 



Lemma 3.2 ([12]) Let < p < q < oo and a function g on I d be represented by the series 

9 = 9k, 9k G Sr(fc), 



/or which 



1/9 

£ 112(^)1*1^,11 | < 



T/ien g £ L q and there holds the inequality 



( 



oo. 



1/9 



|| 2 (i/jp-V«)l*|i 



9k\\ q p 



with some constant C depending at most on p, d. 
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Lemma 3.3 Let < p, 9,q < oo and ip : 
following. 



Then for every f € , there hold the 



(i) For p > q, 



\\f-RA(f)\\ q « 



sup 2~^ k \ 
fceZ^\A 



9 < min(g, 1), 



l/e* 



[Erf\A{2- #) f , 0>min(g,l), 



where 8* :- 



1 



l/min(g, 1) - 1/0' 
(ii) For p < q < oo, 

\\f-RA(f)\\ q « ||/|| B w 

1 

where q* :- 



sup 2 -,/ ' (fc)+(1/p ~ 1/9)|fc|1 , 
fcez^\A 

E fce zt\A{2^ (fe)+(1/p - 1/9)|fc|i r*) 1/9 * 



> q, 



l/q-1/0- 
(iii) For p < q = oo 



#A(/)lloc « ||/|| B W^ 



P. « 



sup 2 -^W+l fc l 1 /P, 
fcez^\A 



A 



< 1 
> 1, 



where 9' :- 



1 - 1/9' 



Proof. 

Case (i): p > q. For an arbitrary / 6 B^g , by the representation (|2.27j) and (|2.19p we have 

ii/--ra(/)iij « E n»(/)ii5 

fceZ|\A 

with any r < min(g, 1). Therefore, if # < min(g, 1), then by Theorem 12.21 and the inequality 

< ||gjfe(/)|| P we get 

1/9 



( 



\\f-RA(f)\\ q « 



E n»(/)ii 
/ 

< sup 2^^ 



1/6 



k€Zl\A 



< ll/lloW SUp 



E {2* W ll«b(/)U' 

\^fceZ^\A 



(3.2) 
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If 6 > min(q, 1), then 

Wf-RA(fW q « E = E V m h k {f)\\ q Y{2-^)y, 

where v = min((7, 1). Since v/0 + = 1, by Holder's inequality with exponents 9/u,0* /u, the 
inequality ||<2%(/)||g < ||<?fc(/)||p and Theorem 12.21 we obtain 



1/8 



RA(f)\\ q « 



E {2^\\q k (f)\\ g } 6 



i/e* 



i/e* 



(3.3) 



This and (|3.2p prove Case (i). 



Case (ii): p < q < oo. For an arbitrary / £ , by the representation (|2.27|) and Lemma 13.21 we 
have 

II. 

Therefore, if 6 < q, then 

\\f-RA(f)\\ q « 



Rami « E {2 (1/ ^ 1/9)|fc|i n%(/)ii P r- 

k&±\A 



E {2 (1/p ~ 1/9) " £|i n%(/)ii P } 6 



1/6 



1/0 



« sup 2 -^(fc)+(l/^-l/^)|fc| 1 J- {2^\\q k (f)\\ p } 6 



k€Z<[\A 



fceZj\A 



If 6 > q, then 



« ||/|| B w sup 2 -^W+( 1 /p-i/'?)|fe|i 



\\f-RA(fw q « E {2 (1/p - 1/9)|fc|i ik fc (/)ii P r 

fcez|\A 

= ^ {2^ (k) \\q k (f)\\p}^2-^ k)+( - 1/p - 1/qWl } q . 
keZf\A 



Hence, similarly to (13. 3D . we get 



( 

| 2 -^(fe)+(l/p-l/ 9 )|fe|i-j.g* 

^fcez^\A 



1/9* 
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This completes the proof of Case (ii). 

Case (iii): p < q = oo. Case (iii) can be proven analogously to Case (ii) by using the inequality 

||/-i?A(/)||oo « £ 2l fe l^|| % (/)|| p . 



Denote by the set of triples (p,9,q) such that < p, 9,q < oo. According to Lemma [331 
depending on the relationship between p,9,q for (p,6,q) £ R+, the error ||/ — i?A(/)||g of the 

approximation of / S has an upper bound of two different forms: either 

\\f-RA(f)\\ q « ||/|| B w sup 2 ~^ + ^- 1 '^\ k \\ (3.4) 

P' e k& d + \A 

or for some < r < oo, 

/ \ iA 



RA(f)\\ q « 



Let us decompose into two sets A and 5 with An B = &s follows. A triple (p, 9,q) € Rj_ 
belongs to ^4 if and only if for (p, 9, q) there holds (|3.4p . and belongs to .B if and only if for (p, 0, g) 
there holds (j3.5|) . By Lemma f3.31 ^4 consists of all (p,0,q) G satisfying one of the following 
conditions 

• P > q, 9 < min(g, 1); 

• P < q, 9 < q; 

• p < q = oo, 9 < 1, 

and B consists of all (p, 9, q) 6 Ri. satisfying one of the following conditions 

• p > q, 9 > mm(q, 1); 

• p < q, 9 > q; 

• p < q = oo, 9 > 1. 

We construct special sets A(£) parametrized by £ > 0, for the recovery of functions / € U^f 
by -Ra(£)(/)- Let < p, 9, q < oo and a € R+, /3 G R be given. We fix a number e so that 

< e < min(a — (1/p — l/q) + , 
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and define the set A(£) for £ > by 
A(£):= < 



'{k £ U{ 
{k G Z d + 



{ a -(i/ p -i/ q ) + )\k\ l+ (3\k\oc <£}, (p,fi,?)eA 

(a -(l/ p -l/ g ) ++E /d)|A;| 1 + ( / g- Olfcloo <£}, (p,e,q)eB, (3 > 0, 
(a -(l/p-l/?) + -E)|A;|i + 09 H-eJIfcloo < £}, (p,9,q)eB, (3<0. 



Theorem 3.1 Ze£ < p, 9, q < oo and a € R+, (3 £R, /3 ^ 0, such that 

l/p < min(a, a + (3) < max(a, a + 0) < r. 
Then there holds true the following upper bound 

sup \\f-R H0 {f)\\ q « 2-«. 



feU 



a,/3 



(3.6) 



Proof. If (p,9,q) € A, by Lemma 13.31 we have 

sup ||/-i2A(0(/)ll« « SU P 2 



-(a-(l/p-l/ 9 )+)|fc| 1 +/3|fc| 00 ) 



We next consider the case (p, 6*, g) G B. In this case, by Lemma |3T3| we have 
sup || 



feu, 



R A(®\\ T q « E 2-^-( 1 /f-V?)+)|fc|i+/9|fc|oo) ) 

fcez*\A(0 



for t = 9*,q*,9'. For simplicity we prove the case (p,9,q) € -B for r = 1 and p > q, the general 
case can be proven similarly. In this particular case, we get 



sup 



f^U 



II/--Ra«)II 9 « E 2- a l fc l^^ fe l- =:E(£). 



(3.7) 



feezi\A(0 



We first assume that f3 > 0. It is easy to verify that for every £ > 0, 



s(0 



w(0 



(3.8) 



where M(x) := maXj- g M Xj for x € M d , and 

W(0 := {i£lf:(a| e/d)(l, x) + {(3 - e)M(x) > £}. 

We put 

V(£, s) := {x e W(0 : £ + s - 1 < a(l, x) + /3M(x) < £ + s}, s € N. 
then from (13.81) we have 



(3.9) 
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Let us estimate |V(£, s)|. Put V*(£, s) := V(£,s) — x*, where x* := {vd)^ l ^\ and v := 
a + /3/d. For every y = x — x* G V*(£, s), from the equation (l,x*) = £/v and the inequality 
a(l, x) + f}M(x) < £ + s we get 

a(l,y) +/3M(y) < s. (3.10) 

On the other hand, for every x G s), from the inequality a(l,x) + /3M(x) < £ + s and 
(a + e/d)(l, x) + (/3 — e)M(x) > £ we get M(x) — (l,x)/d < e -1 s. This inequality together with 
the inequality a(l, x)+j3M{x) > s — 1 gives (1, x) > + ((1 — e~ l j3)s + l)/v for every x E V(£, s). 
Hence, for every y = x — x* G V*(£, s), 

(l,y) > ((1 - ff- 1 ^ + (3.11) 

This means that V*(£, s) C V(s) for every £ > 0, where V'(s) C M. d is the set of all y G M rf given 
by the conditions (|3.10p and (13. llf) . Since V'(s) is a bounded polyhedron and consequently, 

\V(H,s)\ = \V*&s)\ < \V'(s)\ x s d , 

combining (|3.7p and (|3.9|) . we obtain 

oo 

sup \\f-R A (Q\\ q « 2-^2- s S d x 2-*. 



If /3 < 0, similarly to (|3.7p and (|3.8p . we have for every £ > 0, 

E(0 x / 2-( al '*)-^( a! >dz 1 

where 

W'(0 := {x G : (a - e)(l, x) + {fi + e)M(x) > f }. 

From the last relation, similarly to the proof for the case /3 > 0, we prove (|3.6p for the case 
P < 0. □ 

We construct special sets A'(£) parametrized by £ > 0, for the recovery of functions / G U^q 
by Let < p, g, # < oo and a G be given. In what follows, we assume the following 

restriction on the smoothness a of B a Q : 

< a x < a 2 < ... < a d . (3.12) 

We fix a number e so that 

< e < a,2 — a\, 

and define the set A'(£) for £ > 0, by 



where a' 




l/g)+|fc|i < e>, (p,0,?)eA, 
-l/?)+|fc|i <£}, (p,e,q)eB, 
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Theorem 3.2 Let < p, 9, q < oo and a G satisfying the condition (|3.12[) and 

l/p < ai < a d < r. 
Then there holds true the following upper bound 

sup ||/ - R A >(o(f)\\ q « 2^. (3.13) 

Proof. Let us first consider the case (p,0,q) E A. In this case, by Lemma 13.31 we have 
sup ||/--R A '(o(/)ll* « sup 2 -^-^ p - 1 '^\ k ^ « 2-<. 

We next treat the case (p,6,q) G B. In this case, by Lemma |3.3| we have 

sup \\f-R A '(®(fW q « E 2-^ fe )-( 1 / P -i/«) + l^| 1 ), 
J p.» fcez*\A'(£) 

for r = 6*,q*,6'. For simplicity we prove the case (p,9,q) G -B for r = 1 and (l/p — l/g)+) = 0, 
the general case can be proven similarly. In this particular case, we get 



sup \\f - R A , i0 (f)\\ q « J2 =:£(£)■ (3.14) 

£U p> e feezi\A'(e) 



It is easy to verify that for every £ > 0, 

E(0 x / 2-^ x Ux, (3.15) 
./w({) 

where 

W(0 := {i£K{: (a',x) > £}. 

We put 

V(£ s) := {x E W(f ) : £ + s - 1 < (a, x) < £ + s}, s G N. 
then from (|3.15p we have 

oo 

S(0 ^ 2-C^)2-|^,*)|. (3-16) 



Let us estimate |V(£, Put V*(£, s) := V(£, s) — x* , where x* := (ai) £e . For every 

y = x — x* G V*(£, s), from the equation (a,x*) = £ and the inequality (a,x) < £ + s we get 
(a, y) < s and therefore, 

< s/ttj, J e [<*]• (3.17) 

On the other hand, for every x G V(£, s), from the inequality (a,x) < £ + s and (a, x) — e(l', x) = 
(a',x) > £ we get (l',x) < e _1 s, where 1' := (0, 1, 1, 1) G This inequality together with the 
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inequality a\X\ + ad(l',x) > (a, x) > £ + s — 1 gives x\ > £,/a\ + ((1 — e ad)s + l)/ai for every 
x G V(£, s). Hence, for every y = x — x* G V*(£, s), 

2/1 > ((l-e-^Js + lVoi, % - > 0, j = 2,...,d. (3.18) 

This means that V*(£, s) C V'(s) for every £ > 0, where V'(s) C M rf is the box of all y G W 1 given 
by the conditions (|3.17p and t|3. 18j) . Since 

\V(Z,s)\ = \V*fos)\ < \V'(s)\ x s d , 

by (gUD and (f37T6|) . we obtain 

oo 

sup ||/-i? A «)|| 9 « 2^V2' s s d x 2"*. 



4 Sparsity and optimality 

Lemma 4.1 Lei < p, 9, q < oo and a G M+ ; /3 G R, /3 ^ 0, swc/i t/iai 

1/p < min(a, a + /3) < max(a, a + 0) < r. (4.1) 
T/ien there holds true the following asymptotic order 

|G(A(0)| x Yl 2 ' fc|1 x ( 4 - 2 ) 
*eA(£) 

where 

^ = (a + P/d-(l/p-l/q)+, 0>O, 
\a + ^-(l/p-l/g)+, /3<0. 

Proof. The first asymptotic equivalence in (|4.2|) follows from the definitions. Let us prove the 
second one. For simplicity we prove it for the case where p > q, the general case can be proven 
similarly. 

Let us first consider the case (p, 9,q) G B, (3 > 0. It is easy to verify that for every £ > 0, 

V 2l fc l 1 x / 2^dx, (4.4) 



where 



:= {x G : (a + e/d){l,x) + (J3 - e)M(x) < £} 
and M(x) := max Jg ui Xj for x G M d . We put 

V(£, s) := {x G W{0 : £/v + s - 1 < (1, x) < tjv + s}, s G 
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From the inequalities /3 > e and M(x) — (l,x)/d > 0, x E Ri, one can verify that for every 
x S W(£), (l,x) < £/v. Hence, we have 



/ 2^dx < 2*/" V 2~ s |V(£,s)|. (4.5) 



Let us estimate |V(£,s)|. Put V*(£,s) := V(£,s) - x*, where x* := (z/d) From the 

equation (l,x*) = £/v, we get for every y = x — x* € V*(£, s), 

s - 1 < (1,3/) < s. (4.6) 

and 

(a + e/d)(l,y) + (/9-e)M(y) < 0. (4.7) 

This means that V*(£, s) C V'(s) for every £ > 0, where V'(s) C M rf is the set of all y € M. d given 
by the conditions (|4.6p and (|4.7p . Notice that V(s) is a bounded polyhedron and |V'(s)| X ,s d_1 . 
Hence, by the inequality 

\V&8)\ = \V(£,8)\ < \V'(s)\, 
(|4.4p and (|4.5p . we prove the upper bound in f|4.2j) : 

oo 

^ 2 |fc|i <<; 2 C/"^2-s d - 1 X &l v . 

To prove the lower bound for this case, we take k* := [£/dv\l £Zf. It is easy to check k* £ A(£) 
and consequently, 

2 |fe|1 > 2 |fc * 11 » 2 ?/iy . 

fc G A(5) 

The case (p, g) £ B, {3 < can be proven similarly with a slight modification. To prove the 
case (p, 9,q) £ A it is enough to put e = in the proof of the case (p, 8, q) € B. □ 

Lemma 4.2 Let < p,9,q < oo and a € satisfying the condition (13.12P and 

1/p < ai < ad < r. 

Then there holds true the following asymptotic order 

|G(A'(0)| - X] 2 ' fc|1 ~ 2«/( ai -( 1 / p - 1 /' ? )+). (4.8) 

feGA'(0 



Proof. The first asymptotic equivalence in (14.8j) follows from the definitions. Let us prove the 
second one. For simplicity we prove it for the case where p > q, the general case can be proven 
similarly. 
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Let us first consider the case (p, 8, q) E B. It is easy to verify that for every £ > 0, 

V 2l fc l 1 x / 2^ x Ux, (4.9) 



fceA'(f) 

where 



W(0 := {x£K|: (o',x)<a 

We put 

y(£, s) := {x E W(0 : £/a x + s - 1 < (1, x) < £/ai + s}, s E 
One can verify that for every x E W(£), (l,x) < £/oi- Hence, we have 



/ 2( 1 ' x )di « 2 ? / ai V 2- s |y(e,s 



(4.10) 



Let us estimate |V(£, s)|. Put V*(£,s) := V(£, s) - x*, where x* := (ai) 1 £e 1 . From the 
equation (1, x*) = £/di, we get for every y = x — x* £ V*(£, s), 

s-1 < (l,y) < s. (4.11) 

and 

(a',y) < 0. (4.12) 

This means that V*(^, s) C V'(s) for every £ > 0, where V'(s) C M. d is the set of all y E M rf given by 
the conditions (|4.11|) and (|4.12|) . Notice that V'(s) is a bounded polyhedron and |V'(s)| x 
Hence, by the inequality 

mMi = < \v(s)\, 

(|4.9p and (|4.10p . we obtain the upper bound in 



fceA'(f) s=o 

To prove the lower bound, we take k* := L^/aiJe 1 E Z£. It is easy to check k* E A'(£) and 
consequently, 

2 l*li > 2 |fe * 11 > 2 €/ai . 

fceA'(£) 

□ 

Remark The grids of sample points G(A(£)) and G(A'(£)) are sparse and have much less elements 
than the standard dyadic full grids which give the same recovery error. For instance, if we take 
the standard dyadic full grids G(A*(£)) where A*(£) := {k E : ^j/cloo < £} and the number v 
is as in (j4.3jl . it is easy to verify that the linear sampling algorithm i?A*(£) on the grids G(A*(£)) 
gives the worst case error 

sup ||/-i? A *(£)(/)l| g ~ 2_e " 
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The number of sample points in G(A*(£)) is |G(A*(£))| x 2 d ^ u . Whereas, due to Theorem 13.11 
and Lemma 14.11 we can get the same error by the linear sampling algorithm Ra(£) on the grids 
G(A(£)) with the number of sample points |G(A(£))| - ■ 

The following two theorems show that the linear sampling sampling algorithms Ra{£) on sparse 
grids G(A(£)), and Ra'(£) on sparse grids G(A'(£)) are asymptotically optimal in the sense of the 
quantities r n and g n . 

Theorem 4.1 Let < p, 9, q < oo and a G R + , @ G R, /3 ^ 0, swc/i i/iai 

1/p < min(a, a + /3) < max(a, a + /3) < r. 

Assume that for a given n G Z+, £ n zs t/ie largest nonnegative number such that 

|G(A(£ n ))| < n. (4.13) 

TTien R<\(£ n ) defines an asymptotically optimal linear sampling algorithm for r n := fn{Upo)q an d 
Qn ■= Qn(Upf) q by 

RA(t n )(f) = L n (X:,$* n ,f) = Yl E W^Wkj, (4.14) 

fceA(£„) je/ d (fe) 

where X* := G(A(£ n )), <£* := {V'&jIfceAten), jei d (k)> an d there hold true the following asymptotic 
orders 

( n -a-l3/d+(l/p-l/q)+ p > q 

sup^ HZ-^aJH, * r n x ft, x n _ Q _, +(1/p _ 1/3) ' ; (4.15) 

Proo/. 

Upper bounds. Due to Lemma 14.11 we have 

n x 2*"/" x |G(A(£ n ))| < n, 
where ^ is as in (|4.3p . Hence, we find 



n -<x+P/d-(l/p-l/q)+ /3 > 0, 
n -a+/3-(l/p-l/q)+^ /3 < 0. 



(4.16) 



By Lemma l3.1l and ()4.13j) . -Ra(£„) is a linear sampling algorithm of the form (|l.ip as in ()4.14p and 
consequently, from Theorem 13.61 we get 



Qn < r n < sup ||/-flA({„)(/)|| g < 2 in . 
These relations together with (|4.16j) proves the upper bounds of (|4.15p . 
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Lower bounds. We need the following auxiliary result. If W C L q , then we have 

r n {W) q » inf sup ll/U,. (4.17) 

X n ={xi}y =1 Cl d f € w. f{xi)=0, j=l,-,n 

For the proof of this inequality see [T71 Proposition 19]. Since ||/||, > ||/||p for p > q, it is sufficient 
to prove the lower bound for the case p < q. Fix a number r' = 2 m with integer m so that 
max(a, a + 0) < min(r', r' — 1 + 

We first treat the case /3 > 0. Put fc* = fc* (77) := 77I for integer rj > m. Consider the boxes 
J(s) C I d 

J(s) := {x£l d : 2-^ +m Sj < Xj < 2^ +m ( Sj + 1), j G [d]}, s G Z( V ), 

where 

Z(tj) := { S CZ|:0< Sj - < 2^ m — 1, j 6 [d]}. 
For a given n, we find 77 satisfying the relations 

n X 2l fe *l 1 X 2^- m ) = |Z(r/)| > 2n. (4.18) 

Let X n = {x J }" =1 be an arbitrary subset of n points in F\ Since J(s) n J(s') = for s 7^ s', and 
|-^(??)| > 2n, there is Z*{r\) C -£(77) such that |Z*(r/)| > n and 

x„n{u seZ , {v) J(s)} = 0. (4.19) 
Consider the function g* G EJ?(fc*) defined by 

5 * — \2-"\it*h-P\k*\c°+\k 9 h/p ^ M k * jS+r , /2 , (4.20) 

sgZ*(r?) 

where M k * >a + r //2 are B-splines of order r' . Since we have 

||<7*||, x / \2~ a l fc *l 1_/3 l fc *l°° + ( 1 / p_1 /^l A: *l 1 , (4-21) 

and 

\\g*\\ p x A 2- a l fc *l 1 - /3 l fe *l-. 

Hence, by Corollary 12.11 there is A > independent of r\ and n such that g* G U^f . Notice that 
Mk* !S +m-i(x), x G" J(s), for every s G Z*(jj), and consequently, by (|4.19j) g*{x 3 ) = 0, j = 1, ...,n. 
From the inequality (|4.17p f)4.21j) and (|4.18p we obtain 

Qn » \\9*\U x n-*-^ +1 / P -i/«. 

This proves the lower bound of (|4.15p for the case f3 > 0. 

We now consider the case j3 < 0. We will use some notations which coincide with those in the 
proof of the case (3 > 0. Put k* = k*(rj) := (j],m, ...,m) for integer 77 > m. Consider the boxes 
J(s) C I d 

J(s) := {xel d : 2^ +m Sl < x x < 2- r > +m {s l + 1)}, s G Z{rj), 
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where 

Z(rj) := {sGZf :0<si< 2^~ m - 1, Sj = 0, j = 2, d}. 
For a given n, we find r\ satisfying the relations 

n x 2 fe * x 2^"™ = \Z(rj)\ > 2n. (4.22) 

Let X n = {x j, }™ =1 be an arbitrary subset of n points in I d . Since J(s) n J(s') = for s ^ s', and 
1^(^)1 > 2n, there is Z*{r\) C ^(r?) such that |Z*(r/)| > n and 

^n{U sgZ , w J(s)} = 0. (4.23) 

Consider the function g* £ Hf(k*) defined by 

g * : = A2 -(a+/3-l/p)fe| £ M fcv+r , /2 , (4.24) 
sez*(jj) 

where Mj,* iS+r '/2 are B-splines of order r'. Since we have 

\\g*\\ q x A2-( a+/3 - 1 / p+1 /'?) fc i*, (4.25) 

and 

|| 5 *|| p x A2-( Q+ W. 

Hence, by Corollary 12.11 there is A > independent of r\ and n such that g* € Notice that 

Mk* !S +m-i{z), x $ J( s ), f° r every s € Z*(rj), and consequently, by (|4.23p g*(x^) = 0, j = 1, ...,n. 
From the inequality (jlTTl) (|4~25l) and (14221) we obtain 

fti(D3f)* » lb* II, * n-^^VP-i/,. 
This proves the lower bound of (|4.15p for the case f3 < 0. □ 

Theorem 4.2 Let < p,9,q < oo and a £ satisfying the condition (|3.12[) and and 

1/p < ai < ad < r. 

Assume that for a given n £ Z+, £ n is i/ie largest nonnegative number such that 

|G(A'(£ n ))| < n. (4.26) 
Taen -Ra(£„) defines an asymptotically optimal linear sampling algorithm for r n := fn(Up'g) q and 

Qn ■= Qn(Upf)q by 

fceA'(£„) je/ d (fc) 

where X* := G(A'(£ n )), <I>* := {V'ftjlfteA'ten), jei d (k)> an d there holds true the following asymptotic 
order 

SUp ||/-i2 A (C n )(/)ll9 ^ r « ^ ^ „-a 1+ (l/p-l/,) + _ (428) 
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Proof. 

Upper bounds. For a given n G Z + (large enough) , due to Lemma 14.21 we can define £ = £ n as the 
largest nonnegative number such that 

n x 2^ ai -( 1 / p - 1 /< ? )+) x |G( A' < n. (4.29) 

Hence, we find 

2 -e« x n -oi+(i/p-i/ ff ) +> ( 4>30 ) 

By Lemma f3.1l and (I4.26P i?A'(£„) is a linear sampling algorithm of the form (II. ip and consequently, 
from Theorem 13. 131 we get 

ft, < r n < sup \\f - R A i(( n )(f)\\ q < 2"^. 
These relations together with ()4.30p proves the upper bounds for (|4.28p . 

Lower bounds. As in the proof of Theorem 14.11 it is sufficient to prove the lower bound for the case 
p < q. Fix a number r = 2 m with integer m so that < min(r, r — 1 + 1/p). In the next steps, 
the proof is similar to the proof of the lower bound for the case /3 < in Theorem 14.11 Indeed, we 
can repeat almost all the details in it with replacing a + ft by a\. □ 

Remark Concerning the asymptotically optimal sparse grids of sampling points G(A(£ n )) and 
G(A'(£ n )) for r n (Upf) q , Q n (Upf) q and r n (U^ g ) q , Q n (U^ g ) q , it is worth to notice the following. 
Let set A and B be the sets of triples (p,6,q) introduced in Sectional 

(i) For every triple (p, 9, q) G A or (p, 9,q) G B in the case (3 > 0, the grids G(A(£ n )) and G(A'(£ n )) 
were employed in [H [9] for sampling recovery of periodic functions with bounded mixed smooth- 
ness. 

(ii) For every triple (p, 9, q) G A, there is only one family of asymptotically optimal sparse grids 
G(A(£ n )) and G(A'(£ n )). Whereas, for a triple (p,9,q) G B, there are many families of asymptot- 
ically optimal sparse grids G(A(£ n )) and G(A'(£ n )) depending on parameter e > 0, for r n (Upf) q , 
Qn(Vpf)q and r n(Up fi ) q , Q n (U* g ) q , respectively. 



5 Sampling recovery in space B2q 

In this section, we extend the results on sampling recovery in space L q of functions from B^' g in 
Sections [3] and [H to sampling recovery in space Bg tT . 

Lemma 5.1 Let < p,9,q,r < oo, d/q < 7 < min(r, r — 1 + 1/p) and ip : — > M + . Then for 
every f G -8^^ , i/iere /joW the following 



^a(/)|| b?t « 11/11 w ^ 



' sup 2~^ k)+ \ 9<t, 
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Proof. From (I2.2ip there holds true the following inequality for every g of the form (12.20|) . 

\\g\\ q « 2^- 1 /^\ k ^\\g\\ p . (5.1) 

Hence, for every / G B^ e \ by Theorem 12.31 and the inequality 

y - RA(fw BlT « £ {^-'^HQkinWpV . 

By use of this inequality and additionally Theorem 12.31 in a way similar to the proof of Lemma 
I3.3f i) we can prove the lemma. □ 



Let < p, 9, q, r < oo and 0,76 R+, /3 £ 1 be given. We fix a number e so that 

< e < min(a - (1/p - l/q) + , | 7 - (3\), 
and define the set A (£) for £ > by 



A"(0 := 



{k G Z + 
{A; G 

^{A; G Z + 



( a _ (i/p _ l/gJ+JIfcli + ( 7 - ^IfcU < £}, < r, 

(a - (1/p - l/g) + + s/d)\k\x + ( 7 - /3 - e)^^ < £}, > r, /? > 0, 

(a - (1/p - 1/(7)+ - e)|A;|i + ( 7 - /3 + e)|fc|oo < £}, > r, /? < 0. 



The following theorems and lemma are counterparts of the corresponding results on sampling 
recovery in space L q of functions from B^g in Sections [3] and HI They can be proven in a similar 
way with slight modifications. 

Theorem 5.1 Let < p, 9,q, r < oo, a, 7 G R+ and /3 G M, /3 7^ 7, satisfy the conditions 

1/p < min(a, a + 0) < max(a, a + 0) < r, d/q < 7 < min(r, r — 1 + 1/p). 
Then there holds true the following upper bound 

sup \\f-R^ (£ Jf)\\ B , « 2-«. (5.2) 



p,S 



Lemma 5.2 Lei < p, 9, q, r < oo ; a, 7 G R+ and /3 G R, f3 7^ 7, satisfy the conditions 

1/p < min(a, a + 0) < max(a, a + 0) < r, d/q < 7 < min(r, r — 1 + 1/p). 
Then there holds true the following asymptotic order 

\G(A"(0)\ ~ £ 2 ' fc|1 x 

k€A"(£) 

w/iere 

'a-fr-P)/d-(l/p-l/q)+, /3> 7 , 
a-(7-/3)-(l/p-l/g)+, /3< 7 - 
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Theorem 5.2 Let < p,0,q,T < oo, a, 7 £ R_|_ and /5 6 1, /3 / 7, satisfy the conditions 

1/p < min(a, a + /3) < max(a, a + /3) < r, d/q < < min(r, r — 1 + 
Assume that for a given n € Z+, ^ n zs £/ie largest nonnegative number such that 

|G(A" (£„))! < n. 

Then R/\"(£ n \ defines an asymptotically optimal linear sampling algorithm for r n := r n(Up'g)^ T 
and g n := Q n (Upf) B i r by 

Ra" «„)(/) = L n (K,KJ) = E E /(2- fc i)^j, 

fceA"(? n )ie/ d (fc) 

where X* := G(A (£n)), := {V'fc,j}fc g A''(£ n ) je/ d (fc)' anc ^ ^ ere ^ ree the following asymptotic 
orders 



SUP ||/-# A "(£ n )(/)b?. T ~ ^ X 



-a-(/9- 7 )/d+(l/p-l/ g )+ i /3 > 7; 



p,8 



6 Optimal cubature 

Let X n = {x J }™ =1 be a set of n sample points in I d , A n = {Aj}™ =1 a sequence of n numbers. For 
a / E C(I d ), we want to approximately compute the integral 



/(/):= / f{x)dx 
Jl d 

by the cubature formula 

n 

I n (X n ,A n ,f) := E A i/(^)- 
3=1 

For C C(I d ), to study the optimality of cubature formulas for / £ VF, we use the quantity 

:= inf sup \I{f)-I n (X n ,A n ,f)\. 

X n ,A„ f£W 

Every linear sampling algorithm L n (X n , $ n , •) of the form (|l.ip generates the cubature formula 
I n (X n ,A n ,f) where 

A n = {Aj }"■=!, Aj = / (fj(x) dx. 

Jl d 

Hence, it is easy to see that 

\I(f)-I n (X n ,A n ,f)\ < \\f-L n (X n ,$ n J)\\i, 
and consequently, from the definitions we have the following inequality 

in{W) < r n {W)L (6.1) 
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Theorem 6.1 Let < p, 9 < oo and a £ R+, /3 G R suc/i that 

1/p < min(a, a + (3) < max(a, a + f3) < r. 
Assume that for a given n 6 r L + , £ n is the largest nonnegative number such that 

\G(A{£ n ))\ < n. 

Then i?A(g n ) defines an asymptotically optimal cubature formula for iniUpf) by 

I n (X* n ,A* n ,f):= Yl X kJ(2- k j), 
2-*jex* 

where 

X*:= G(A(£ n )), := {h,jh-k je x* > x k,j ■= J ^k,j{x)dx, 
and there hold true the following asymptotic orders 

( n -a-P/d+(l/p-l)+ R > 

W)-W,K,f)\ - in(U«f) x | n _ a _, +(1/p _ 1)+) ' j; <a (6.2) 

Proof. The upper bound of (|6.2[) follows from (|6.ip and Theorem 14,11 
To prove the lower bound of (|6.2h we observe that 

in(W) > inf sup |/(/)|, 

^n={^}! j ' =1 Cl <1 f eW: f( x i)=0, j=l,...,n 

and for the functions g* given in (|4.20p and (|4.24|) we have 

itf) = Ibli- 

Hence, we can see that the lower bound is derived from the proof of the lower bound of Theorem 
HU □ 

In a similar way, we can prove the following 

Theorem 6.2 Let < p,9 < oo and a £ satisfying the condition (|3.12p and ai > 
Assume that for a given n 6 Z + , £ n is £/ie largest nonnegative number such that 

|G(A'(e„))| < n. 

Then i?A'(£„) defines an asymptotically optimal cubature formula for in(Up 6 ) by 



I n (X*,A* n J):= Yl X ^f( T 



k 3) 
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where 



X* := G(A' (£„)), A* := {X k j} 2 -k jeX *, X k j = J ifj k ,j(x)dx, 
and there holds true the following asymptotic order 

sup \I(f)-I n {A* n ,X* n ,f)\ x i n (U a pfi ) x „-«!+(!/?-!)+. 

Remark If in Theorems 16. II and 16.21 we assume 1 < p < oo, then 

i (IT*) - [n-*-^, /?>0, 
n[ p > 8 ' ~ \n—P, f3<0, 

and 

i n (U» e ) x n- a K 
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